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Abstract 



Q ■ Here is presented a four-body potential model for q'^cf systems 

Qs . which includes both the spin and flavour degrees of freedom, extend- 

"^ I ing the formalism presented already in the spin independent situation. 

Q_i' This allows an application to a realistic situation, which is chosen 

to be KK scattering. It is seen that because of the gluonic effects 
D I in this multi-quark system, the KK attraction resulting from the 

i-^ ' quark-exchange mechanism gets appreciably decreased compared to 

that emerging through the naive two-body potential approach. 
V ; PCAS number(s):12.40.Qq, 13.75.Lb, 12.38.Lg, 14.40.Cs 



1 Introduction 

With strong evidence in favour of quark confinement both from experiment 
(failure to find free quarks) and theory [|l], |^, we have some understanding of 
the quark-quark interaction for large distances. Using this, along with the un- 
derstanding of the short distance quark-quark interaction obtained through 
*E-mail address (Internet) :MASUD@PHCU.HELSINKI.FI 



perturbative QCD, different models of the quark-quark interaction have been 
used such as the MIT bag model |P, ^ and the constituent quark poten- 
tial model. In the constituent quark potential model the quark- ant iquark 
interaction is represented by a potential which is well motivated by QCD 
for small distances (given by the one-gluon exchange mechanism), but is 
modified so as to incorporate the confining potential as the limit for large 
distances. Moreover, the current masses of the QCD lagrangian are replaced 
by effective masses, termed constituent masses, which are fitted, along with 
other parameters of the model, to experimentally known quantities related to 
some set of hadrons. For quarks of known (fitted) masses interacting through 
a space dependent potential, one can set up and solve a Schrodinger equation 
for their dynamics. In this way the constituent quark potential model, when 
improved to incorporate relativistic effects, explains in a consistent and uni- 
fied way most of the observed mesonic states, from the pion to the upsilon, 
as quark-antiquark states with different values of orbital and radial quantum 
numbers 0. 

But that is not enough. A successful model of strong interactions should 
be able to describe also possible systems having three or more quarks/antiquarks. 
It is not clear yet how, if at all, the quark potential model can be applied to 
multi-quark systems (with more than three quarks). Perhaps the simplest 
approach is to take the many-body hamiltonian as a sum of hamiltonians 



corresponding to all pairs of particles involved, the basic method used in 
d) 0) H) P • But this has many theoretical as well as phenomenological (such 
as the van der Waals force problem [|l^) flaws. Keeping this in mind, a 
four-body potential model for a quark-exchange mechanism in g^g^ systems 



was proposed in |Tl| in a spin independent situation, taking into account the 
effects of the gluonic degrees of freedom as well in a non-trivial way. Only 
for small distances does this agree with the sum of two-body (Fj ■ Fj) po- 
tentials model. For, roughly speaking, interquark distances greater than 0.5 
fm it qualitatively agrees with the flux tube model [O, O of the adiabatic 



surfaces of the gluonic field. 



As shown in |jT^ , there are three flux tube topologies for the Ng = Ng = 2 
system, each of which will determine the ground state in a different region 
of configuration space. Since these are linearly independent of each other, 
we need a basis containing at least three gluonic states to describe the glu- 
onic field of this multiquark system. This is different to what the two-body 
potential model says for the g^g^ system; the two-body model would need 
only a two dimensional colour basis |^. Thus in [|rT[] the two-body potential 



model was written in a redundant colour basis 

|l)e = |li3l24)e, |2), = II14I23), and |3)^ = iSiaSsi),, (1.1) 

corresponding to the three basic states in the flux tube model, and then the 



suggested changes, forming the proposed model, were made. The (confine- 
ment) potential part of the hamiltonian of the system was written as 



^(91929394) 



E F. ■ F,v. 



ij, 



with 



i<j 

Crl + C. 



:i.2) 

:i.3) 



This quadratic, rather than the theoretically better motivated coulomb-plus- 
linear, form of the two-body potential was used basically for computational 
convenience. Similarly for the kinetic energy part, the non-relativistic expres- 
sion was used, also for simplicity and numerical convenience. These features 
have been retained in the present work as well. 

The overlap matrix A^ in this redundant basis is as given by eq.(2.8) of 
[pj]l , the potential matrix by eq.(3.10) there, and the matrix element of the 
non-relativistic kinetic energy operator between any two states is expressed 
through eq.(3.11) of the same (these are also given by eqs.( [L.4|) , (|1.6|) and 



LSD below respectively, with / = 1). The model proposed in [|r^ gives the 



following expressions for these matrices: 

/ 1 



N ^ N{f) 



-f 
1 



If 
J 



:i.4) 



,{X'\H\X), - 

,{X'\K\X), = N{f)]l?^(j: 



g{X'\K\X)g+ g{X'\V\X)g, 



1/2 



2r7i,- 



mm: 



x 



and 



:i.5) 



V 



-\{Vi^ + V2i) 
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f 12 + VM 
-^13 - ^24 
— f 14 — f 23 

-|K4 + ^^23) 



2/ 
3v^ 



/ -2(t;i3 + ^24) \ \ 



2/ 



symmetric 



\ 



+^^14 + ^^23 

-f 12 - f 34 

2{Vi4 + V23) 

-f 24 - f 13 
2(t^l2 + ^^34) 
3 I +^^13 + ^24 
+Vr4 + V2^ 

+ \Df{l-f) 



J 



1.6) 



The basis now is actually |l)g, \2)g and |3)g, instead of |l)c, |2)c and |3)c; 
the new subscript g refers to the gluonic degree of freedom, instead of the 
colour degree of freedom represented by the subscript c. H is the hamiltonian 
of the system, and N{f)x',x is the gluonic states overlap factor g{X'\X)g. 
The matrices in the above three equations are achieved by multiplying the 
two-body potential model expressions for the off-diagonal elements of overlap, 
kinetic energy and potential energy matrices respectively by a factor / which 
tends to 1 for all rij <^ &^^ and to when any Tjj becomes ^ K^ ■ ^^ 
the proposed model in |n|, the diagonal elements were the same as obtained 
in the two-body model calculations, apart from the |-D/(1 — /) term in the 
3,3 element of the proposed potential matrix (eq. (|1.6|) ). The reasons for the 
introduction of this term are described in the discussion from eq.(4.18) till 



the end of section 4 of flTTl . It should be clear from the discussion there that 



the model is actually defined in a basis which is not redundant even for all 



Tij vanishing. 

This treatment of the diagonal and the off-diagonal elements of the two- 



body model based matrices was motivated by the work presented in |1J, [15 
along with the similarity for large interquark distances in space dependence 
of the diagonal elements of the K, V and A^ matrices in both the two-body 



potential and the flux tube model. The calculations reported in [T^, |15 
actually show that in the flux tube model the coupling of the three gluonic 
states |1) , |2) and |3) decreases exponentially with inter-quark distances. 
For the space dependent factor / multiplying the off-diagonal elements, the 



choice used in [^ was 



/ = exp(-^^4) , k = -kh, (1.7) 



with k a numerical coefficient. This is the simplest choice from a compu- 



tational point of view. An alternative choice of / is suggested in |]16|, [1 



Both of these forms have been studied in [|18| , |T9|, ^, which aim at extract- 
ing the gluon field overlap factor / from a calculation using lattice Monte 
Carlo techniques. This work is in progress, but so far neither of these two 
candidates is preferred by the calculations. 

Developing a formalism, based on the model proposed above, describing 
the dynamics of meson-meson interactions also requires taking into account 
the "slower motion", the quark position dependence of the wave function. 

6 



This was done, for the spin-independent case, in sections 5 and 6 of the 
previous work [|ll|. There we used our model to specify a theoretically refined 



(but non-relativistic) hamiltonian of the q^(f system in our gluonic basis, 
and solved approximately for the wave function using the 'resonating group 
method' common in nuclear [^ and also recently in particle physics [^ 



The approximation corresponded to specifying parts of the quark position 
dependent wave function before solving the Schrodinger equation for the 
rest. 

In the present work we introduce the spin and flavour degrees of freedom, 
and as a first application apply the formalism to a physical meson-meson 
system, namely KK. This is done in section 2. The solution for the total 
wave function of the system thus obtained gives us numerical results for the 
corresponding meson-meson phase shifts for the elastic as well as the non- 
elastic meson-meson scattering, along with a condition for the existence of 
a bound state of the whole system. These results are reported in section 3. 
This is followed by our conclusions. 

2 The model applied to KK systems 

In this section we present the formalism for a realistic situation. This means 



that along with the gluonic basis described in [11|, we have to deal with 



the spin basis as well. Moreover, flavour dependence has to be considered. 
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The hamiltonian of the system, written now in this spin-gluonic basis, would 
also include the hyperfine term. To complete the formalism, we have to 
incorporate also the quark position dependence of the wave function. With 
quark contents of the q^(f system like that of KK , the flavour wavefunction 
can be written generally as Isls, with / standing for a light (up or down) 
quark. We label these four particles 1,3,4 and 2 respectively. Thus the pair 
(2, 3) would be composed of strange quark and antiquark, and (1, 4) of light 
quark and antiquark. In this way each of the particles 2 and 3 has a mass 
higher than of those belonging to the other pair (1 and 4) by a ratio which 
is the same as that of the strange quark mass irig to the up (or down) quark 
mass m. This mass ratio is denoted by s in this paper. Anti-symmetrization 
of the total wave function is not necessary in the present case since we do 
not have any two identical fermions. 

For the spin dependent part of the basis we use, as in the quark model, the 
states arising through the spins of the quarks only. So any pair of particles 
(quark or antiquark) would have a total spin of one or zero, thus forming a 
spin triplet or singlet respectively. This also means that the total spin of the 
whole q^q^ system can have a value of zero, one or two. Being interested in 
the ground state of the J^ = 0+ sector of the system, we focus on the spin 
states with the total spin of the system as zero. This is consistent because, 
as is mentioned later, tensor and spin-orbit forces are neglected in this work. 



Hence, the hamiltonian separately conserves L and S i.e. the total orbital 
angular momentum and the total spin of the system. This allows us to restrict 
our considerations to the S = sector, as mentioned above, meaning, in turn, 
that we will be dealing only with L = spatial wave functions. 

In every one of the three channels of the previous section (correspond- 
ing to the three gluonic states ll)^,, \2)g and |3)c,), our four particles can be 
grouped into two mesonic sub-clusters. Each of these clusters may have a 
combined spin of zero or one and hence the q^q^ system may be composed of 
either two spin singlets or two triplets. This means that there may be two in- 
dependent spin channels for each of the three gluonic channels above. Thus, 
there are six independent states of the system in hand. The corresponding 
six spin states are written in the notation of Appendix D of (see Appendix 
A of the present paper for details) as: 
In the first channel (with the gluonic part of the base state as |l)g): 

\lS), = \P,^P2i)s and |lT), = |Vi3-V24).. (2.1) 

In the second channel: 

\2S)s = \PuP2l)s and |2T), = |Vi4 ■ V23).. (2.2) 

And in the third channel: 

|35), = |5i2%). and |3T), = |Ai2 ■ A34).. (2.3) 



In this notation Sij and Ajj stand for the scalar and axial vector spin wave 
functions respectively, and the pseudoscalar and vector spin wavefunctions 
Pij and Vjj are defined in terms of their linear combinations. 

Except in the diagonal terms corresponding to the PP (i.e. pseudoscalar- 
pseudoscalar) sector of the second gluonic channel to be discussed below, the 
flavour wave functions are taken to be trivial everywhere, just giving rise 
to an isospin conserving factor as the overlap of any two of them. ; This 
happens in the absence of any mechanism for flavour change of a quark or 
antiquark. Actually, flavour changing is possible in any channel through 
annihilation of quarks and antiquarks of the same flavour. Our consideration 
of these processes in the PP sector of the second channel only is, therefore, 
an approximation. This sector is singled out because it is here that the 
annihilation effects are apparently most significant and cannot be neglected 
in any realistic model of the processes involving mesons. This is because 
they are supposed to be responsible for the mass difference between the 
(pseudoscalar) isoscalar and isovector mesons. The annihilation effects are 
negligible in the (spin) vector-vector sector of the second channel, because 
of the small difference in mass of the spin one isoscalar and isovector mesons 
i.e. u and p. 

When these annihilation processes are incorporated, the fiavour wave 
function gets mixed with the quark position dependent part. This follows 
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because the size of a mesonic cluster depends upon the masses of the quarks 
(antiquarks) it contains. Thus the PP sector of the second channel is con- 
sidered separately below. For all other channels, we can consider the quark 
position dependent part of the wave function-referred to as 'quark wave 
function' in the following-separately from the flavour part. This quark wave 
function is a function of four 3- vectors ri,r2,r3 and rj. These can be re- 
placed by their combinations, with one of them as the overall centre-of-mass 
co-ordinate of the whole system R^ and three others which are taken, here, 
to be different in different channels. Writing explicitly, these are: 
In the first channel (with the gluonic part of the wave function as \V)g) 

R-i = 7-, , yi = ri - rg and Zi = r2 - rj. (2.4) 

1 + s 

In the second channel 

ri + rj — r2 — rg 
R.2 = , y2 = ri - VI and Z2 = r2 - rg. (2.5) 

And in the third channel 

ri 4- r2 — rg — rj 
R3 = — , yg = n - r2 and Z3 = rg - r^. (2.6) 

1 + s 

The quark wave function in any channel is written, in the following, as 
a product of two factors, one being a function of R^ and the other of y^ 
and Zfe only, for k = 1,2 or 3. The former is denoted by XkiiJ^k)-, with / 
designating the spin state (singlet-singlet or triplet-triplet), and the latter 

11 



by C.k{yk)Ck{2k), with ^fc(yfc) and Cfc(zfc) corresponding to the two mesonic 
clusters of the channel k. The spatial dependence of these on y^ and z^ is 
taken to be gaussian in consistency with the choice of the quadratic form of 



the inter-quark potential in eq.( |1.3| ). But the Xfc/'s a^re treated as variational 



functions to be determined by solving the approximate coupled Schrodinger 



equations, using the 'resonating group method' |^ . 

With the above mentioned forms of the quark wave functions in different 
channels, the total state vector of the whole g^g^ system is written as 

l^(?i, 92, ga, Qi, g)) = J2 \k)g\kI)s\k)f'i/jci'Rc)Xkii'Rk)^kiyk)Ckizk), (2.7) 

kl 

with 

^''^y''') = T^—Wrm^M-yl/'^dli] and Cfc(zfc) = -— — -^exp[-z^/442]- 
(2vrrf^i) (2vrrf^2) 

(2.8) 

This is actually the case with annihilation neglected. In that situation the 
mesons represented by ^i,Ci5^3 and ^3 have one light and one strange (anti- 
strange) particle. On the other hand, that denoted by .^2 has both particles 
as the light ones (up and down) and (2 has both the quark and antiquark as 
the heavier ones. It follows thus from the properties of the solutions of a 3-d 
harmonic oscillator that dii,di2,d^i and ^32 have a particular value, say, d', 
d2i has a different one, say, d, and ^22 differs from all these having a value 
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denoted by d" . Quantitatively 



d'^ lm{m, + m) Is + l ^ d"^ [2^ /T ,^ ^, 

It must be pointed out that we are using the approximation of neglecting 
the spin dependence of the size of any cluster. For the absolute magnitudes 
of the sizes, the relation d"^ = v3-R^/2 is used relating the radius of a meson 
composed of the light mesons only to the r.m.s. charge radius -R„=0.6 fm 
of a nucleon whose qqq wave function is generated by the same quadratic 
confining potential. 

The PP (or S) sector of the second channel needs special consideration 
because of the annihilation and creation processes, making the pairs (2, 3) 
and (1,4) mixtures of ss,uu and dd in the flavour space. The particular 
combinations depend upon the physical mesons taking part in the scattering 
process which may be rj and/or rj'. In the isovector sector only (2, 3) has an 
amplitude for going to, say, uu from its original flavour state as ss, resulting in 
1] or 1]'. As the size, and hence the quark wavefunction, of a particular cluster 
is related to the masses of the quark it contains, these amplitudes for having 
different flavour contents imply that in this particular channel we cannot 
write the quark and flavour wavefunctions separately as it is possible for the 
other channels. Rather these two are mixed here, giving us the combined 
quark-flavour wavefunctions (except X2s(R-2)) as 

13 



|25)/, = |Mi4)/,|M23)/, 



with 



V'li) Ms) fg 


for 


v'v' 


mesons 


Vu)fq ^23) /9 


or 






\v[i)fq\V2-3)fq 


for 


T]r]' 


mesons 


IVU) fq\V23) fg 


for 


riT] 


mesons 


I'^li) fq V23) fq 


for 


Tcrj' 


mesons 


. I'^li) fq V23) fq 


for 


TCT] 


mesons 



(2.10) 



Here 



Ur^. 



\Vij)fq 

Wij)fq 
kli)fq 



COS 6 



sin 6 



\dd + uu) f ^ , , . /ii _v ^ / N 
^6(r»i) - sine\ss)fC2{rij) 



_V2 
\dd + uu) j 



V2 
vr)/6(ri4) = k)/6(y2). 



6(rij) +cos6'|ss)/C2(ri 



(2.11) 
(2.12) 



1 

72 n3/4 



1 



/o ^2 N^/- exp [-4/4^2,] and C2(r.,) = - 
(27rc?^i) (27rrf22 



'2 A 3/4 



exp [-47442], 
(2.13) 



where ^21 = '^,'^22 = d" , and 9{= 34.7°) is related to the mixing angle 
Op = —20° of flavour singlet and octet resulting in i] and rj' (see pages III. 68- 
69 of [H). 

After this discussion of the total wavefunction, we write the Schrodinger 
equation for the system in hand: 



{H - E,)\^{qi,q2,q3,q4;g)) = 



(2.14) 



14 



where H is the total hamiltonian and E^ is the total centre-of-inass energy 
of the q^(f system. This also means that the overlap of {H — Ec)\^) with an 
arbitrary variation |(5^) of the state vector |^) vanishes. In \5'^) we consider, 
as in resonating group method calculations, only the variations in Xki- Thus 
we write 

kiw' 

M'Rc)6xkii'Rk)Uyk)Ckizk) f{kU{ki\g{k\H-E,\i)g\ij)s\i)fM'Rc)xiA'RiMyiKii^i) = O- 

(2.15) 
To do the four space integrations implied in the overlap of the diquark di- 
antiquark position dependent {H — -Ec)!^) and |5\l/), any of the three sets 
of three- vectors defined by eqs.( |2.4| ),(p.5|) and ( p.6| ), along with Re, can be 



used. The choice Rfc,yfc and z^ has, however, a clear advantage. The arbi- 
trary variations 6xki(Rk)'s for different (but continuous!) values of R^, and 
of k and / as well, are linearly independent and hence their coefficients in 
eq.( |2.15| ) should be zero. With the trivial Re integration performed to give 



a finite result using, say, box normalization, this leads to 



^ / d'ykd'zkUyk)Ck{zk)f{k\s{kI\s{k\H-E,\l)g\lJ),\l)fXiA'RiMyiM^i) = 0, 

(2.16) 
for /c, /=1,2 or 3 and I, J = S or T (except the special case of k, 1=2 and 
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I, J = S where a mixed quark- flavour wavefunction is used). 

In ||Tl| (and also in the previous section) the suggested hamiltonian of the 
system was mentioned in the spin independent context. Despite the addi- 
tional degrees of freedom incorporated here we proceed in the same fashion, 
considering first the hamiltonian in the two-body potential model limit and 
then modifying the off-diagonal elements. Now the spin dependent part of 
the hamiltonian, composed of the terms corresponding to hyperfine (contact 
as well as tensor) and spin-orbit interactions, has to be considered as well, 
along with that representing the annihilation effects. But our job here is 
simpler than this because of our constraint to the S-wave ground states, im- 
plying that the only additional spin-dependent term in the hamiltonian to 
be dealt with is that representing the hyperfine contact interaction. 

For the hyperfine term in the two-body potential model limit we take the 
expression given by one gluon exchange (used, with some modifications, in 
[^]), and sum over all the pairs: 

yhyp = y: V^y^ = - E F^ ■ F,|^^53(r.,) S, • S,. (2.17) 

Numerical values of a*-', fitted to the light meson spectroscopy below, will be 
taken as varying with the sum of the masses of the particles i and j; thus 
each of them will eventually be replaced by a", a'-f or a'^/, I standing for a 
quark(antiquark) of light mass (i.e. up or down) and s for a strange (heavier) 
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one. 



The pair annihilation and creation effects (where considered) are repre- 
sented by a hamiltonian term V"' (denoted by Ha in the hamiltonian appear- 
ing in ) operating in the ff avour space only. This should be a sum of two 
terms belonging to the pairs 14 and 23. Thus 

V^« = 1/^1 + 1/25 • (2.18) 



The matrix elements of V/^ are written as 



f{uu\Vr^\uu)j = f{dd\Vr^\uu)f= f{dd\Vr^\dd)f = l, (2.19) 

f{uu\V^j\ss)f = f{dd\V^^^\ss)f = Vl^ and f{ss\V^^^\ss)f = n. (2.20) 

in the corresponding flavour spaces. The mass (or flavour) dependence shown 
here is in qualitative agreement with the mass dependence of the annihila- 
tion term Ha appearing in [Q. In the above, / and n are phenomenological 
parameters to be fitted to the masses and flavour wave functions of vr, 77 and 
7]' mesons. 

This specification of the hamiltonian of the system means that H in 
eq.(P3q) is to be replaced by fsT + 1/^ + 1/<^ + eLi ^h with V^ = V^f + V'^yp 



i.e. a sum of the confinement and hyperfine potentials. Using eqs. (|1.2|) and 
( 2.17 ) for V^^ and V^^^ respectively, it can be seen easily that now the matrix 



elements of V^ in our spin basis, in the two-body potential model, would be 
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given by 



,{X\VP\X'), 
{Vij)x„x'^ 



i<j 

V,, ,{X\X')., ''™° 



with 



(2.21) 



53(r,,).(X|S,-S,|X%.(2.22) 



UillUj 



In this form, s{X\Vp\X')s is very similar to the expression for V appearing 
in eq.( |1.2|) . So, as far as its matrix elements between the gluonic states 
are concerned, the whole formalism developed in ||1T| (also reported in the 
previous section) can be utilized, provided the spin state dependence of the 
matrix coefficients Vij, replacing Vij, of Fj ■ Fj is taken care of. Thus, the 
matrix element of the V^ term in the hamiltonian between any two gluonic- 
spin states appearing in eq.( |2.1(j| ) is given by (see eq.(|1.6|) for comparison) 



/ 



yp 



-1(^13 



Vo 



24^1, 1 



9 



^12 + Vm 

— Vi3 — V24 

—Vii — V23 



2/ 



1,2 



-1(^14 + V^23)2,2 ih 



symmetric 



-2(V^13 + V^24) \ \ 

+VI4 + Vrs 

/ 2(^4 + ^3) 
+Vu + Vu 

V -V^24-V^13 

2(Fi2 + Vm, 
+Vr3 + V21 
+Vu + ^3 



2,3 



V 



+IJD/(1-/) 



3,3 

(2.23) 



in the gluonic-spin basis ll)^!!^)^, \l)g\lT)s, ■ ■ ■ , |3)c,|3T)s. 

The above matrix should give 6 x 6 = 36 matrix elements of the V^ 
operator in its 6-dimensional basis. But as written above, it has only nine 



elements. Actually, every term in every matrix element of the above is meant 
to stand for a 2 x 2 matrix defined by 

( {yij)kS,lS iVij)kS,lT \ 



{Vij)k,, 



(2.24) 



\ iVij)kT,lS (yij)kT,lT J ' 

for k, 1=1,2 or 3. (yij)kLu is given through eq.( p.22| ) and / multiplying the 



D term in 3,3 element of eq.( p.23| ) is a 2 x 2 identity matrix. To determine 
the overlap of any two spin states and the corresponding matrix element of 
the operator Sj ■ Sj, use has to be made of the definitions given through 
eqs.(p.l| ) ,( |2.2| ) and ( |2.3| ) above. The results are reported in Appendix A. 

As far as the spin dependence is concerned, the other terms in the hamil- 
tonian are unit operators. Using the results reported in the previous section 
for the matrix elements between gluonic states, we write eq.( p.l6| ) as 



ij •' 



;^/ 



hyp 



K,ki,ij(Rk, H-0 + Vfcj ;j(Rfc, Rz) + V^fuCRk, Hi) 



-\E,-J2^i] A/'fc/,zj(Rfc, R/ 



i=l 



xu(Ri 



0,(2.25) 



for I = 1,2 and 3 along with J = S and T. This gives six equations, each for 
one of the six possible values of the pair of variables k, I with k = 1, 2 or 3 
and I = S or T. Here lCki,Ui V^{ u, V^^u ^"^^ ■^ki,ij are defined through the 
following equations: 

y'd3R;/Cfc,,;j(Rfc,R;)xzj(R;) = y'c^VfcO?=^Zfeefc(yfc)G(zfc)i^fc/,axzj(R.oeKy/)0(zO 

(2.26) 
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(2.27) 

(2.28) 

J d^R[Afki,iA^k, R'i)xij{^'i) = J d'ykd^ZkUyk)Ck{zk)Nki,ijXij{^iMyi)(iizi] 

(2.29) 
with Kki^ij, V^j ij and Vj^j^fj representing the matrix elements of the K, V^^ 
and V^"^^ operators between the spin and gluonic states appearing in eq.( p.l(: ). 



Nkidj is the overlap of these states, calculated using the results mentioned 
in the previous section along with those in Appendix A for the spin overlap 
factor s{kl\lj)s. 

The spatial integrations appearing on the RHS of the eqs. (|2.26| ) to ( p.29| ) 



are done after substituting the expressions for Kki,ij, V^j u, Vf^j^fj and Nki^ij, 
obtained through the procedure outlined above. In these calculations the 
diagonal {k = I) and off-diagonal {k ^ I) cases were dealt with separately. In 
the former case, xui^i) is linearly independent of the integration variables 
Yk and z^ and thus was simply taken out of the integrations. In the case of 
off-diagonal terms (with k ^ I), the integration variables (y^ and z^) were 
replaced by their equivalent combinations with one identical to R^, and the 
other one independent of it. Integrating out the vector independent of R^, 
expressing the remaining one on the RHS (apart from R;) in terms of R^ 
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and R/, we got the results for ICki,ij, V^j u, Vf^j^u and N'ki,u after comparing 
with the LHS of the corresponding equation. 

Where annihilation is considered, we have to use the combined quark and 
flavour wave function \2S)fq defined through eq.( p.lO|) . Thus in the diagonal 



term corresponding to the 25* channel we have, in place of the k = I = 2 and 
I = J = S term in eq. (|2.16|) . 



d'rud'r^s fg{2S\,{2S\g{2\H ~ E,\2)g\2S)s\2S)fg X2s(R2) 



= Jdh,id''rr3 f,{2S\ s{2S\g{2\K + VP- {E, - ^m,)|2)g|25),|25)/,X2s(R2) 

1=1 

+ Jdhudh,-,f,{2S\VTS)fgX2s{n2). (2.30) 



As expressed through eq.( p.lO|) , the form of 125)/^ depends upon the physical 
content of the 2S channel. This would result in different expressions for 
each of IC2s,2S,^2S2S^^2s'^2S ^^^ J^2S,2S for different pairs of mesons in the 
channel. In the following calculations we restrict ourselves to just the lowest, 
in threshold energy, channels: rji] in the isoscalar, plus rjir and 77V in the the 
isovector sector. This is done because of our special interest in the behaviour 
of the KK system near the threshold (see the next section). 

The results thus obtained for these non-local kernels, for all the values of 



/c, /, / and J, appearing in eq.( p.25|) are reported in Appendix B. Substitution 
of these in eq.( |2.25D would gives six coupled equations. However we neglect 
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all the connections to the third gluonic channel, justified to some extent by 
the absence of any significant effect of removing the third channel in the 



spinless case (see fig.4 of [|TT1). This leaves us with just four equations (two 



of these are written below as eg. ( p. 331) and ( |2.34]) ). The off diagonal (i.e. with 



k y^ I) terms in these equations tend to zero for large inter-cluster distances. 
Thus for consistency with the observed meson spectroscopy we would require 
the constant term in each of the diagonal parts to be identical to the sum 
of masses of the mesons present in the corresponding channel. Fitting, in 
this way, to the masses of K, r], r]', vr, K*, u (or p) and mesons, we get the 
following values of the above mentioned free parameters of the formalism: 

m = 277 MeV, s = 1.955, C = 456 MeV, 
a^^ = 1.583, a[' = 1.561, af = 1.501, 

/ = 272 MeV and n = 67 A MeV. (2.31) 

For C (see eq. (|1.3|) ), the equality of kinetic and potential energies of a har- 
monic oscillator is used, giving us 

C = — ^-4 = ^ = -270 MeV/fm^. (2.32) 

After this parameter fit (except for k to be discussed in the next sec- 
tion), we write down two of the four coupled equations mentioned above. 
The remaining two equations would involve vector mesons. These are not 
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incorporated beyond this stage because of our above mentioned neglect of 
channels opening at energies significantly higher than the KK threshold. 
1 



Mk + Mk 



+eo / d^K2 



1 1 

2m 6 



V^. - E, 



(giiR? + quKl + gio) + C(Ri) 



exp 



— e2R.2 — ciR,]^ 



G'(Ri,R2)iJ X25(R2)=0, 



and 



(2.33) 



Ma + Mb- V^^ - E, 



X2S(R2 



+eo / rf^Ri 



1 1 



~^l (921R? + g22R2 + g2o) + C{Ki) 



exp 



-e2R2 — ciRi 



-G'(Ri,R2)i/ Xi5(Ri)=0, 



(2.34) 



for KK <-^ ah, where a and h are the two mesons in the second channel 
{r]ri,riTT or i]'7r). C(Ri) and G'(Ri,R2) appearing in the above are 



C(Ri) 
G(Ri, R2) 



i (feiR? + &o) - ^ {e, + lc- 2m{s + 1) 



(2.35) 



/loexp -(ei + e[ - lii)Rl - (e2 + /i2)R^ 



X 



a^^ exp (/13R1 ■ R2) + ctg exp (-/13R1 ■ R 
+/2oexp -(ei + e'l + /2i)Ri - e2R2 

a" exp (/22R?) + sttf exp (-/22Ri^ 
+/3oexp [-(ei + e[ + hi)Rl - (e2 + /32)R2 x 

a[' exp (/33R1 • R2) + a'f exp (-/33R1 ■ R.2)l • (2.36) 
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Moreover, 

6 3m2s(27rrf2)3/2- ^ ' > 

It should be noted that in the coefficients of V^^ and V^^' ^^^ reduced 
masses of the two pseudoscalar mesons of the particular channel now appear. 
This is done so as to ensure that the terms involving V^^ or V^^ Si'^^ ^^^ ^'^^~ 
rect kinetic energy of the relative motion of the interacting physical mesons. 
Other symbols appearing in the above equations are defined in Appendix B 
at appropriate places. 

The kernels of the integrals appearing in the off-diagonal parts in the 
above coupled equations contain non-separable parts 

exp(/i3Ri • R2) ■ ■ ■ exp(/33Ri ■ R2). The presence of these "non-separable po- 
tential" terms makes the solution of the coupled equations rather involved. 
To avoid that complication, we can solve our problem in the approximation 
of replacing these terms by their truncated expansions which would leave 
us with an inexact but easily manageable form of the equations. With that 
strategy in mind, the above equations were solved ffist for the case of no 
hyperfine interaction (thus avoiding non-separable terms) by setting H = 
in the above two coupled equations. The method used for that is explained 
below for the full interaction case. The resulting phase shifts for this no 
hyperfine case (reported in the next section) are so small that it would be a 
good approximation to take the variational wave functions in the absence of 
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hyperfine interaction as the wave-functions corresponding to a freely propa- 
gating plane wave. Using this approximation for the variational wave func- 
tions xislRi) and xislR-i) even in the presence of hyperfine interaction, we 
looked for a reasonable separable approximation to our non-separable poten- 
tial terms. As far as the terms exp(/i3Ri ■ R2) and exp(— Z13R1 ■ R2) are 
concerned, it was seen to be a very good approximation to just replace them 
with the exponential expansion up to the second power in /13R1 ■ R2. But 
the terms multiplying /30 are not so easy to manage. Their separable ap- 
proximation involved two parameters which had to be adjusted for each of 
value of the kinetic energy in the centre of mass frame. Written explicitly, 
our approximation has been to take 



^30 exp -(ei + e'l + /3i)R? - (e2 + ^32)R2 



X 



[a[' exp (-/33R.1 ■ R2) + a[' exp (/33R1 ■ R2)} ^ 
nih^Q^^ |exp -Ti (ei + e'^ + /31 + 62 + h2 - hz) (R-? + R.2 

- -ri (ei + e'l + /31 + 62 + ^32 + ^33) (R? + R2)] } > (2-38) 

in (^(Ri, R2), and hence in the above two coupled equations. Here rii and ri 
are the above mentioned two energy dependent parameters. Actually these 
are to be used in eq.( |2.33| ) only. Those appearing in eq.( p.34D are different 
and, therefore, are denoted by n2 and T2 instead. 
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Both sides of eg. ( P^B]) , multiplied by Rlx2s{R2) and integrated over R2, 



were plotted as functions of Ri. This comparison showed that in this way 
even the worst discrepancy could be reduced to less than 10 percent of the 
total hyperfine coupling for that particular value of -Ri and the on-shell mo- 
mentum J5c(2). After calculating the T and 5* scattering matrices, the results 
were checked for any possible deviation from unitarity of the S matrix and 
symmetry of the T matrix (required by "reciprocity" of inelastic scattering, 
see, for example, p. 528 of [Q). The discrepancy was for some cases as bad 
as 30 percent, implying that the above approximation needs to be improved 
by, for example, iterating it many times. This improvement remains to be 
made, although it can be easily shown that this re-adjusting of the values 
of ni,n2,Ti and T2 (with improved functional dependences of xis'(R'i) and 
X25(R'2)) is not needed for the range of energy where our immediate interest 
lies (e.g. below the KK threshold in the first channel). This follows because 
the formal momentum space solutions (see eqs. (p.39|) and ( p.40| ) below) of 



the above coupled equations, for that range of energy, would be of the form 

xip) = — . , -. X Constant. 

Changing values of ni,Ti,n2 and T2 in that situation would just affect the 
constant coefficients of . ] . and -it-t-^, leaving the momentum, and hence 

Ai(pi) A2(p2)' *= ' 

the space, dependence of the solutions Xis and X25 unchanged. 
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With the above approximation the integrand appearing in the couphng 
terms in both of our coupled equations (|2.33| ) and ( p.34|) are products of 
two factors, each of them is a function of Ri or R2. This means that in 
this form the two coupled equations can be solved exactly, using the method 



demonstrated in appendix B of [|TT]. The procedure used in the present 
case was to first write eqs.( 2.33|) and ( p.34 ), in the approximate form (see 
eq. (|2.38| )), in momentum space. For incoming waves in the first channel, the 
formal momentum space solution of these equations would be 



XisiPi 



X2S{P2 



5(Pl-Pc(l)) 



pia) 



l(l) 



1(1) 



(1) 



QY'A^ie^) + Q'2>B2{e2) + Q's'A^ie^ + lu) + Ql" B^ie^ + I 



Ai(Pi 
+Qi^^A2inei + e'l + /31 + 62 + ^32 - /ss) 

+Qi^^A2{nei + e'l + /31 + 62 + h2 + hs) 



l(2) 



(2) 



grMi(ei)+Qr'5i(e 



)(i) 



12 J 



22yl 



tlly 



A2(P2 

+Qf /li(ei + e'l + /21 - ^22) + Qf ^i(ei + e[ + hi + I 
+Qf''Ai{ei + e[ - In) + Q^ Bi{ei + e[ - I 
+Qf^Ai{T2ei + e'l + /31 + 62 + /32 - /ss) 
+gf Ai (rsei + e'l + /31 + 62 + /32 + hs) 



The new symbols appearing in these equations are defined in Appendix C 



(2.39) 



(2.40) 
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It is to be noted that in the above equations pi and p2 have been re- 
placed everywhere by pi and p2 respectively, utilizing the spherical sym- 
metry of our problem. Multiplying eq. (|2.39| ) by plFa{pi,ei) ,plFb{pi,ei), 



pjFa (pi, ei + e[ + hi - I22) ,plFa {pi, ei + e[ + hi + ^22) ,P?-^a {pi, ei + e[ - hi), 
pjFb {pi, ei + e[ - hi) ^PJFa (pi, raCi + e[ + hi + 62 + h2 - hs) 



and pfFa ipi, Ti^i + e'^ + /31 + 62 + ^32 + ^33) in turn and integrating w.r.t. 
Pi gives us 8 equations (Fa(pi,x) and Fb(pi,x) are the Fourier transforms of 
exp [— x-R^] and _R^exp \—xF3{\ respectively). Similarly multiplying eq.( |2.40| ) 
by the Fourier transforms p\Fa {p2, 62) , pjFb {p2, 62) , pjFa {p2, 62 + /12), 
plFb (p2, 62 + /12) , plFa {p2, Tici + c'l + /31 + 62 + ^32 " ^ss) and 

P^Fa {P2, nei + e'l + /31 + 62 + /32 + ^33) and integrating w.r.t. p2 gives us 
6 more equations. These 14 equations can be written as a matrix equation 

QUi = U2 (2.41) 
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with 



U, 



An 



Fb{pc{l),e,) 
Fa{Pc{l),ei + e[ + l2i-l22) 

Fa{Pc{l),ei+e[+l21+l22] 

i^a(Pc(l),ei + e;-/n) 
Fh{pc{l),ei + e[-hi) 



Fa (Pc(l), 7261 + e'l + hi +62 + 1 



32 — '33 



/3 



Fa (pc(l), T2ei + e'l + /31 + 62 + I32 + hs , 









(2.42) 



/ 



Q a 14 X 14 matrix containing many integrals, and Ui a vector containing 



^(62), ^2(62), ■ ■ ■ , A2in6i + e'l + /31 + 62 + h2 + I33) 



and v4i(ei), -Bi(ei), ■ ■ ■ , v4i(r2ei + e'l + /31 + 62 + ^32 + ^33) as its elements. In- 
verting the matrix Q gives these 14 elements of the Ui vector. With these 
values in hand, all the quantities in the expressions for the variational func- 
tions Xis{pi) and X2s{P2) are known. So these can be now simply obtained 
by making the usual replacement of pi and p2 by their on-shell values Pc(l) 
and Pc(2), defined by eqs. (|C.10D and ( p.ll| ), respectively in eqs.( p.39D and 

(HD. 



From eqs. (|2.39|) and ( p.40|) the two T matrix elements Ti^i and T2^i, pro- 



portional to the coefficients of the non-relativistic Green operators 



1 



Ai(pi) 



and 



1 

'A2(P2) 



respectively, can be read off. These are reported in Appendix 
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D. Similarly, for incoming waves in channel 2, the use of U2 as 



/ 



U, 



An 











































FaiPci2), 


e2) 


Fb{Pc{2), 


62) 


Faipci2),e2 


+ h2 


Fb{pc{2),e2 


+ 112 



\ 



Fa ipc{2), nei + e'l + hi + 62 + I32 - I: 



33, 



(2.43) 



\^ Fa {pc{2), nei + e'l + /31 + 62 + /32 + /33J / 
instead of that given by eq.( p.42| ), gives the two T matrix elements T2^2 and 
Ti 2 reported also in Appendix D. 

For the total energy in the centre of mass frame above the higher thresh- 
old, both of the channels would be open. Thus for incoming waves in either 
of them, there would be a loss of flux in the incoming channel (with the total 
flux remaining conserved). Representing this inelasticity by a factor e^, for 
A; = 1 or 2, we can write 



..2i5i 



^1,1 = 1 - 2zri,i = eie"*"^ 

<S'2,2 = 1 - 2iT2^2 = £26^ 



(2.44) 
(2.45) 



For elastic scattering, ei or £2 would be unity for incoming waves in channel 
1 or 2 respectively. 
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Below the lower threshold the situation would be qualitatively different 
as, with both Pc(l) and Pc(2) being imaginary, 6{pi —pc{l)) and 6{p2 — Pc(2)) 
would not contribute to the integration over all the real values of pi and p2 
performed to arrive at eq.( p.41|) . Thus all the terms collected in the vector 



U2 would be absent, leaving us instead with 

QUi = 0. (2.46) 

A non-trivial solution of this equation for the elements of the vector Ui 
requires 

det Q = 0. (2.47) 

giving us a condition for the existence of a bound state of the whole system. 

3 Results 

As should be clear from the expressions reported in Appendix B, our coupled 
equations can describe a number of meson-meson systems. Amongst these 
systems we have chosen KK, keeping in mind that it has been investigated 
by other groups using different models for the quark-quark interaction. An 
important issue is whether the whole KK system has a bound state just 
below the KK threshold or not. According to our method of answering this 
question, the condition for a non-trivial solution of the Schrodinger equation 
to exist for the total energy of the whole system below the lowest threshold 
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is the one given by eq.( |2.47| ). The the Q matrix mentioned there is actu- 



ally a complicated function of the parameters of the formalism fitted above, 



except /c, and the total energy of the whole system, k (see eq. (|1.7] )) is the 
phenomenological parameter of our model of gluonic effects. Our numerical 
calculations, including calculation of the determinant of the Q matrix, were 
done for three values of k in turn. The value k = corresponds to a two-body 
potential model hamiltonian. On the other hand, k = 1/2 fm~ is emerging 



from the lattice gauge theory calculations |T^ for rectangular configurations 



of quark positions. For other configurations, indications |^ are that the 
spatial decrease of gluonic topologies overlap may be slower, and thus we 
have also used an intermediate value (k = 1/6 fm~ ). Which of these, if any, 
would simulate the "experimental" (lattice gauge theory calculations based) 
behaviour of the gluonic overlap is yet to be seen. 

Our numerical calculations showed that for any value of k, the above 
condition for the existence of a bound state (see eq.( p.47|) ) is not satisfied 



for any value of energy below the KK threshold. This is the situation in 
the isoscalar as well as in the isovector sector, whereas in the latter case 
all connections to the nrj channel are neglected as those would not affect 
the answer to the main question being discussed here. On the other hand, 
using a closely related model Weinstein and Isgur |^ get KK bound states in 
both the isoscalar and the isovector sectors, and conclude that the two scalar 
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meson resonances /o(975) and ao(980) can be explained as loosely bound 
KK states. Their model corresponds, in some approximation (see the first 
paragraph of the next section), to ours in the limit k = 0. Therefore it is 
interesting to see if, in the corresponding limit, we can get their results by 
varying our parameters. Our calculations showed that, in the A; = limit, 
we need to multiply our total couplings of the two channels by a factor of 
2.715 before we can get bound states in both the isoscalar and isovector 
sectors. Alternatively, we can get these bound states by multiplying only 
the hyperfine couplings by a factor of 1.926. We get bound states without 
hyperfine interaction as well, but for that an increase by a factor of 3.06 in 
the remaining couplings, with fc = 0, would be needed. This is one of the 
indications in our work that the hyperfine coupling is the main interaction 
arising through quark exchange, and that the hyperfine and other couplings 
arising through confinement etc. have opposite signs. 

Our modification to the two-body potential model proposed in this paper 
(equivalent to using non-zero values of k) implies a decrease in the KK 
coupling. That means that we need to increase the couplings even more so 
as to get binding. The factors so needed for the different values of k are 
reported in table (|37T|), along with the corresponding energy values for the 
resulting binding. 

In addition we report below (in tables (|3.2|) to ( p.TD ) the KK phase shifts, 
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k 

in fm"^ 


hyperfine coupling 


increase 
in the total 
coupling to 
get binding 


increase 

without 

hyperfine to 

get binding 


increase to 
get binding 


B.E. in the 
isovector 
channel 


B.E. in the 
isoscalar 
channel 





1.926 


<1 MeV 


35 MeV 


2.715 


3.06 


1/6 


3.00 


<1 MeV 


53 MeV 


4.26 


6.57 


1/2 


6.43 


<1 MeV 


87 MeV 


8.89 


19.13 



Table 3.1: Increases in the couplings necessary to get binding. 



in the elastic as well as in the inelastic region, for the hyperfine interaction 
increased by the above mentioned factor of 1.926, for the different values of 
k. With that increase in the coupling we get bound states of the whole KK 
system for fc = in both the isoscalar and the isovector sectors, and what is 
explored here is just the effect of our proposed modification to the four-body 
potential. The numerical procedure to get these phase shifts was based on 
eqs. (p.44| ) and ( p.45|) . Each of them is a complex equation and hence can 
be solved for the two quantities e^ (the inelasticity factor) and 6^ (the phase 
shift), with A; = 1 or 2, for each value of energy. 

Also mentioned are the values of the inelastic phase shifts for the incom- 
ing waves in the other channel i.e. rji] in the isoscalar channel and tti]' in 
the isovector one. Moreover, the values of all of these phase shifts in the ab- 
sence of hyperfine interaction (without any increase of coupling) are reported. 
However, it must be emphasized that because of the various approximations 
which we have used in this work, it would be improper to take these numbers 
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(total cm energy) 


full coupling 


without hyperfine 


'^1 
(degs) 


ei 


^2 

(degs) 


5i 
(degs) 


ei 


^2 

(degs) 


997.0MeV 


102.09 


0.729 




1.11 


1.00 




1042.0MeV 


91.49 


0.719 




3.05 


1.00 




1092.0MeV 


96.56 


0.726 




4.05 


1.00 




1095. 5MeV 


99.01 


0.728 




4.25 


1.00 




1142.0MeV 


91.74 




155.94 


3.94 




1.35 


1192.0MeV 


83.41 




147.26 


3.39 




1.46 



Table 3.2: The isovector sector phase shifts for k = Q. 

as precise results of our model. One of the indications of this inaccuracy is 
the violation of unitarity resulting from our separable approximation to the 
actual non-separable terms in the integrands appearing in our coupled equa- 
tions (|2.33|) and (|2.34|) . As mentioned in the paragraph following eq. (p.38|) , 
this approximation affects badly our results for the phase shifts, although not 
our conclusions regarding KK binding. For the elastic region this unitarity 
violation manifests itself in the reported (see tables ( ^^ to (^3)) deviation 
from unity of the inelasticity factor ei. 



4 Conclusions 

In this paper a formalism has been developed to deal with meson-meson sys- 
tems ihaving their dynamics resulting through quark exchange effects, and 
applied, as a first application to a realistic case, to KK systems. Here we 
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(total cm energy) 


full coupling 


without hyperfine 


6i 
(degs) 


ei 


62 
(degs) 


61 
(degs) 


ei 


62 

(degs) 


997.0MeV 


157.18 


0.953 




2.29 


1.00 




1042.0MeV 


125.92 


0.804 




6.56 


1.00 




1092.0MeV 


118.85 


0.779 




9.61 


1.00 




1097.5MeV 


119.76 


0.781 




10.73 


1.00 




1142.0MeV 


112.10 




142.78 


9.26 




4.38 


1192.0MeV 


105.26 




127.53 


7.33 




4.45 


Table 3.3: The isoscalar sector phase shifts for A; = 0. 


E, 
(total cm energy) 


full coupling 


without hyperfine 


6i 
(degs) 


ei 


62 
(degs) 


61 
(degs) 


ei 


62 
(degs) 


997.0MeV 


5.13 


0.997 




0.21 


1.00 




1042.0MeV 


17.10 


0.972 




0.59 


1.00 




1092.0MeV 


30.15 


0.919 




0.82 


1.00 




1095.5MeV 


33.03 


0.905 




0.86 


1.00 




1142.0MeV 


30.84 




179.46 


0.83 




0.27 


1192.0MeV 


27.19 




177.86 


0.74 




0.30 


Table 3.4: The isovector sector phase shifts for A; = 1/6 fm"^ 


E, 
(total cm energy) 


full coupling 


without hyperfine 


6i 
(degs) 


ei 


62 
(degs) 


61 
(degs) 


ei 


62 
(degs) 


997.0MeV 


17.21 


0.972 




0.40 


1.00 




1042.0MeV 


50.39 


0.821 




1.18 


1.00 




1092.0MeV 


79.08 


0.725 




1.79 


1.00 




1097.5MeV 


87.37 


0.715 




1.98 


1.00 




1142.0MeV 


86.34 




151.52 


1.86 




0.90 


1192.0MeV 


78.69 




140.20 


1.56 




0.96 



Table 3.5: The isoscalar sector phase shifts for k = 1/6 fm 
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(total cm energy) 


full coupling 


without hyperfine 


Si 
(degs) 


ei 


S2 

(degs) 


Si 
(degs) 


ei 


S2 

(degs) 


997.0MeV 


0.58 


1.000 




0.02 


1.00 




1042.0MeV 


1.96 


1.000 




0.06 


1.00 




1092.0MeV 


3.38 


0.999 




0.09 


1.00 




1095.5MeV 


3.62 


0.999 




0.10 


1.00 




1142.0MeV 


4.15 




1.25 


0.10 




0.03 


1192.0MeV 


4.25 




1.51 


0.09 




0.04 



Table 3.6: The isovector sector phase shifts for /c = 1/2 fm 



Ec 
(total cm energy) 


full coupling 


without hyperfine 


'^1 
(degs) 


ei 


S2 

(degs) 


Si 
(degs) 


ei 


S2 

(degs) 


997.0MeV 


1.17 


1.000 




0.04 


1.00 




1042.0MeV 


4.18 


0.998 




0.12 


1.00 




1092.0MeV 


8.01 


0.994 




0.19 


1.00 




1097.5MeV 


9.22 


0.991 




0.22 


1.00 




1142.0MeV 


9.28 




3.83 


0.21 




0.10 


1192.0MeV 


8.36 




4.39 


0.19 




0.12 



Table 3.7: The isoscalar sector phase shifts for k = 1/2 fm' 
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had to increase our resulting coupling by some numerical factor before we 
could get a bound state of the whole system, even in the two-body potential 
limit. The variational calculations based on the two-body potential reported 
in ^ claim to get bound states of the whole KK system, concluding that 
the two scalar meson resonances /o(975) and ao(980) can be explained as 
loosely bound KK states in the isoscalar and the isovector sectors, respec- 
tively. Our detailed model of meson-meson dynamics, even in the two-body 
potential model limit, is different to theirs, mainly because of our restricted 
(i.e. only in the 25" diagonal term) incorporation of annihilation effects. This 
neglect of the annihilation effects may have appreciably decreased the KK 
binding arising through our model. This is expected as the annihilation part 
of the hamiltonian, incorporating the process KK ^^ tttt, was ||2^ mainly 



responsible for the KK binding in the calculations reported in ^. Moreover, 
we might be underestimating the hyperfine interaction by treating this inter- 
action partially as a perturbation, although in our work as well the hyperfine 
coupling turned out to be the main interaction arising through the quark 
exchange mechanism. It is difficult to say more about this problem unless 
a more refined treatment of the hyperfine interaction, along with the anni- 
hilation effects, is carried out. On the other hand, the fitting of the model 
parameters in ^ includes adjusting the ranges and normalization of their 
effective meson-meson potentials in an ad-hoc way, and it is not clear how 
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that affects their results. 

Leaving these issues to some future work, we looked for any possible 
change in one of our parameters so as to get KK bound states in the limit 
where our model would roughly correspond to that used in (i.e. for k = 
0), and then determined, with the changed parameter, the effects of going 
beyond that limit i.e. of using our theoretically improved four-body potential 
rather than the naive two-body potential. This means the use of non-zero 
values of k in our terminology. This investigation showed the same trend 
as observed in the spin independent case reported in fTTI]: increasing k, 



i.e. decreasing the gluonic states overlap, results in a significantly weaker 
meson-meson interaction. This means that if we get a KK bound state in 
the two-body potential model limit, we do not necessarily get one with our 
Q CD- inspired refinement of the g^g^ potential. 

Much improvement in the calculations can be made by going beyond 
the approximations we have used, giving quantitatively more precise results. 
But even without this being carried out, this work clearly indicates that the 
theoretical refinement of the four-body potential results in an appreciable 
decrease in a major part of the meson-meson interaction — enough to cast 
doubt on any result based on a naive two-body potential model. 
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Appendix A: The spin basis 

In the main part of this paper, we use the spin basis given through the 
eqs. (p.l|) , (|2.2| ) and (p.3|) . The notation used in these equations is that of the 



Appendix D of 0. In this notation, an orthonormal spin basis is |S'i25'34)s 
and |Ai2 ■ Agj)^, defined through eqs.(D3) and (D4) there. Our remaining 
spin base states are given in terms of these by eqs.(D5-D8) of the same. 
These equations give easily the overlaps of our spin base states, written as 
s{kl\lj)s in the results mentioned in Appendix B, as elements of the overlap 
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matrix 



A 



y symmetric 
in the basis \lS)s, \'^T), 



/l/i -^3/4 yi/i ^ ^ 

1 -^374 -^1/4 ^3/4 -^1/4 

1 -^ y^ 

1 -^3/4 -^1/4 





1 



(A.l) 



/ 



|3T),. 



For the matrix elements of the Sj ■ Sj operators, for different values of the 
indices i and j, in our spin basis, we also used the results expressed through 
eqs.(D9-Dll) of the Appendix D of 0. Some of the results obtained in this 
way are reported here: 



,{IS\ 



rsi 


S21 


Si 


S3 


Si 


S4 
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S3 


S2 
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S4J 
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24 
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4 
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S3 


Si 
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[IT\ 



Si ■ S2 




[T ' 

V 16 


Si • S3 
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S2 ■ S3 
S2 ■ S4 


ll^).= 


^V 16 
+ VT6 




_ S3 ■ S4 _ 




V 16 
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(A.4) 
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.S3 


S4. 



,{2S\ 



\2S), 



"Si 


S2" 


Si 


S3 


Si 


S4 


S2 


S3 


S2 


S4 


.S3 


S4. 



\2S), 



,{2S\ 



Si ■ S2 




"+r 


Si • S3 




-l 


Si-Sj 
S2 ■ S3 


lS)s = 


:i 


S2 ■ S4 
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(A.5) 



(A.6) 



Appendix B: The kernels of the integro-differential 
equations 



The results for (in general non-local) kernels appearing in eg. ( p .251 ), cal- 
culated using the procedure outlined in the two paragraphs following this 
equation, are: 



J^ki,kj — Sij5(Rk — R-fc)' 



for all values oi k I and J. 



(B.l) 



K,ki,kj — Sij5(Ilk — Rfc, 



-[uJki + ujk2\ - 7;^Vl 



(B.2) 



2m—"-\ 
for all values of k, I and J, except those corresponding to the 25* diagonal 

term. Here 

hk 



^kl 
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2mdl^ 
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/i — /s 
9i = 93 = h = h 

/2 
92 

h2 



s + 1 
s + 1 

s 
s + 1 

2s 



(B.4) 



dki and dk2, for A; = 1, 2 or 3, are the same ones which appear in eq. (|2.8| ). 



ycf 


= 5jj5iR^- 


-R'l) 


-^C-AC[dl, + dl,] 
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-- S oiT. Similarly, (B.5) 


^2I,2J - 


= 6jj6{K2- 


-K) 


-p-4C[dl, + dl,] 




(B.6) 



Vstsj = 5/j5(R3 - R'3) X 
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s-1 
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(See the discussion before eq.( |2.9|) for d' 
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,(A;/|/J),iVO 
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for A; = 1 or 3 (B.14) 
for / = 1 or 3 (B.15) 
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the new definitions used are 
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As mentioned in the text, in the case of the diagonal term corresponding 
to the 25* channel, the expressions depend upon the physical content of this 
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channel. Thus we have: 
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(Here f^ = 1, /^^ = 1(1 + 1/s), /f = I/5 
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51 



^[4i + ^2 sin' + ul^ cos' 6] - ^{f^ sin^ 6 + /^^ cos' ^)V^^ 



[B.77) 



Ws,2si^V') 
^2S$siVV) 

^2slsi^V) 

y2lU^v') 

^2S,2siVV) 
^2S,2si^V') 



5(R2 - R2) 

6{K2 - R'2) 
S{R2 - K) 



Svr 



3m^s 

Stt 
3m^s 

Svr 



<5(R2 - R2 
<5(R2 - R2 
<5(R2 - R2 



— C - 2 X 4:C[dl^ cos' 9 + 4^ sin^ 9] 

~C- AC[dl^ + rf^i cos' 9 + d^2 sin' ^] 
o 

--C - 4C[rf2^ + (i^^ sin^ 9 + ^^2 cos^ 9] 



g(a,^^W + afW)cos^g 
(27rrfi,)3/2 



a 



+ a 



23W)sin2( 



sa 



14(0 ^^23(0 pog2 ^ 



+ 



sa: 



+ 



a 



s(27rrfi2)3/2 

23(s) gin2 Q 



?)m?s 
2 X 5(R2 - R2 



(27rdii)3/2 (27rrfii)3/2 s(27rrf22)3/2 



sa 



14(0 



+ 



+ 



(27rrfii)'/' (2vrrfii)'/' s(27rrfi2)3/' 
2/ cos^ 9 - 2 v^2/r2 cos 6^ sin 6* + ra sin^ 9 



5(R2 - R2) 



2/ cos^ 6^ - 2 V 2/ra cos 9sm9 + n sin^ 6* 



5(R2 - R2) [2/ sin' ^ + 2^2/^ cos ^ sin ^ + n cos' ^ 



Appendix C: The definitions used in the mo- 
mentum space solutions 



(B.78) 
(B.79) 
(B.80) 

(B.81) 

(B.82) 

(B.83) 
(B.84) 
(B.85) 
(B.86) 



The following definitions are used in writing the momentum space solutions 
( |2.39|) and (|2.40| ) of the coupled equations. 



Q\ 



(1) 



1 gii 

2m 6 



2 



Fb{Pi,ei) + 



'2m 6 2 



K 
6 



Fa{Pi,ei] 
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Q? 



Q? 
Qf 
Qf 



-Ha%QFa{j>i, ei + e[ + hi - I22) - HsafkoFaipi, ei + e[ + ki + I22) 

^qi2 

2m 6 

-2HaiHwF,{pi,ei + e[-h 



-Fa{pi,ei) nonumber 



ni. 



--Ha'^koll^Fbipi, ei + e[ - In] 



-Hriia^'hoFaipi, nei + e[ + /31 + 62 + h2 - hs) 
Hma^fkoFaipi, TiCi + e'l + /31 + 63 + I32 + ^33) 



1 ?22 j;, / N 



2m 6 2 6 



FaiP2,e2) 



'2m 6 2 



Fa{P2,e2] 



-Ha%oFa{j)2,e2) 
-Hsafl2QFa{p2,e2) 
-2Hai'hoFa{p2,e2 + h2) 



--Ha%olfsFb{p2, 62 + /12) 



-Hn2a'l3oFaiP2, ^"261 + e'l + Z31 + 62 + ^32 - ^33) 
Hn2alH3oFa{p2, T2ei + e'l + /31 + 62 + ^32 + ^33)- 



(C.l) 



(C.2) 



(C.3) 



K = E, + -C -2m{s + l) 
o 



Ak{x) 
Bk(x) 



Co / (i'^Rfcexp 
Co / (i'^Rfcexp 



XJrv?, 



XJrvr, 



Xfcs(R-fe) 



(C.4) 
(C.5) 
(C.6) 
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Xks[Pk) 



d Rfe 



exp [ipk ■ Rfc] XhsCR-k 



(C.7) 



(27r)3/2 

for k = 1,2. Fa{pk,x) and Ff,{pk,x) are similar FTs of exp [— xR^] and 
R^ exp [— xR^] respectively. Moreover 



Ai(pi) 

A2(P2) 
Pc(l) 

Pc(2) 



Pi 



KK 



2/X 
2/^afe 



+ Mk + Mk - E^ - ie 



+ Ma + Mb-Ec-ie 



2fiKK{E, -Mk- Mk) 
2ij,h{Ec-M,-Mh). 



(C.8) 

(C.9) 
(CIO) 
(C.ll) 



Appendix D: The elements of the T matrix 



Consistent with our definition of the T matrix (see eqs. (|2.44|) and (|2.45|) ), the 
four elements of the 2x2 T matrix are (these can be read off from eqs .( 12.391) 
and (g^): 



vr 



i(i) 



i(i) 



i(i) 



i(i) 



Ti,i = 2/i^^-p,(l) [QY>A,{e2) + Q'2>B2{e2) + Q^ A2{e2 + lu) + QT B2{e2 + /12) 



+Q\ M2(riei + e'l + /31 + 62 + ^32 - ^33 
+g^'^A2(riei + e'l + /31 + 62 + /32 + l 



33 j 



TT ,.. IV2 ^ (2) . , . , ^(2) 



T2,i = 2/^,,-p,(l)J-^[gr'Ai(ei) + gr^5i(ei) 

+Qf ^i(ei + e'l + /21 - ^22) + gr^i(ei + e^ + /21 + ^22) 
+Q(')Ai(ei + e'l - /n) + Qf B^{e^ + e'^ - /n) 



1(2) 



+Q7 ^1(^261 + e'l + /3I + 62 + ^32 - / 



33 j 



(D.l) 
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l(2) 



+QI Mi(r2ei + e'l + /31 + 62 + h2 + I 



33 J 



(D.2) 



^2,2 = 2/i,,^p,(2)[gf)/li(ei) + gf5i(ei) 



1(2) 



+QfAi{ei + e'l + hi - I22) + QT Ai{ei + e'^ + /21 + /: 



(2) 



22 J 



1(2) 



l(2) 



+QrMi(ei + e'l - /n) + g[i'^5i(ei + e'^ - hi] 

+Q?Ai{T2ei + e'l + /3I + 62 + ^32 - /3 



f33j 



)(2) 



+Q^'Mi(r2ei + e'l + /31 + 62 + ^32 + ( 



Ti: 



33 j 



^ '- /^i foW^ ^. ^ , r^mn (^ \ , ^W /i ^. , / ^ , riW 



(D.3) 



2/ii^x77Pc(l)W- Qi ^2(62) + Q^^52(e2) + gr'^2(e2 + /12) + Qt' B2{e2 + /12) 
^ V f 2 '- 



1(1) 



+Q^'M2(riei + e'l + /31 + 62 + ^32 - ^33) 
+QSi'^A2(riei + e'l + /31 + 62 + /32 + ^33) 



with pi and p2 in Qi ■ ■ • Qg (see eg s . (10.21) and ( p.3|) ) replaced by vSX) ^'^'^ 
Pc(2) respectively. It was checked numerically that Ti^2 = ^2,1, fulfilling the 
requirement of "reciprocity" in an inelastic scattering (see p. 528 of p^). 
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